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ABSTRACT 

Context. It has been shown that convection in red supergiant stars (RSG) gives rise to large granules causing surface inhomogeneities 
together with shock waves in the photosphere. The resulting motion of the photocenter (on time scales ranging from months to years) 
could possibly have adverse effects on the parallax determination with Gaia. 

Aims. We explore the impact of the granulation on the photocentric and photometric variability. We quantify these effects in order to 
better characterize the error possibly altering the parallax. 

Methods. We use 3D radiative-hydrodynamics (RHD) simulations of convection with C05B0LD and the post-processing radiative 
transfer code 0PTIM3D to compute intensity maps and spectra in the Gaia G band [325 - 1030 nm]. 

Results. We provide astrometric and photometric predictions from 3D simulations of RSGs that are used to evaluate the possible 
degradation of the astrometric parameters of evolved stars derived by Gaia. We show in particular from RHD simulations that a 
supergiant like Betelgeuse exhibits a photocentric noise characterised by a standard deviation of the order of 0.1 AU. The number of 
bright giant and supergiant stars whose Gaia parallaxes will be altered by the photocentric noise ranges from a few tens to several 
thousandths, depending on the poorly known relation between the size of the convective cells and the atmospheric pressure scale 
height of supergiants, and to a lower extent, on the adopted prescription for galactic extinction. In the worst situation, the degradation 
of the astrometric fit due to the presence of this photocentric noise will be noticeable up to about 5 kpc for the brightest supergiants. 
Moreover, parallaxes of Betelgeuse-like supergiants are affected by a error of the order of a few percents. We also show that the 
photocentric noise, as predicted by the 3D simulation, does account for a substantial part of the supplementary 'cosmic noise' that 
affects Hipparcos measurements of Betelgeuse and Antares. 

Key words, stars: atmospheres - stars: supergiants - astrometry - parallaxes - hydrodynamics - stars: individual: Betelgeuse - 

1 ■ Introduction previous work by iLudwi^ (l2006h has shown that effects due to 

. t r , ^ . . . n u^^j the granulation in red giant stars are not likely to be important 

The main goal of t he Gaia mission d Perryman et al. | |200lt ^ ^^^^^^^ -^^^^ 

i LindegrenetalJ M) is to determine high-precision astromet- j^^^ supergiant (RSG) stars are late-type stars with masses be- 
ne parameters (i.e., positions, parallaxes, and proper motions) ^^^^^ ^ j^^^^ ^^^^^^ temperature rang- 
for one bilhon objects with apparent magnitudes in the range ^^^^ 3450 ^^5^ ^^qq ^ ^^^^^ luminosit ies in the range 
Z u- . multi-band and multi- 200O to 300000 Lo, and radii up to 1500 Ro (iLevesque et all 
epoch photometric and spectrocopic data will allow to recon- ^ ^^^^^ luminosities place them among the brightest stars, 
struct the formation history, structure and evolution of the ^j^^^j^ 1 distances. Based on detailed radia tion- 
Galaxy. Among al the objects that will be observed, late-type hyd rodynamics (RHD) simulat ions of RSGs (iFrevtag et al.l2002l 
stars present granulation-related variability that is considered, in iFrevtag & HSfned l2008h. IChiavassa et alJ (l20OTn Paper I 
this context, as noise that must be quantified in order to better hereafter) and ICMavasietaLl (l20ici (Paper II hereafter) show 
characterize any resulting error on the parallax determination. A ^^^^ ^j^^^^ ^^^^^ characterized by vigorous convection which 
Send offprint requests to: A. Chiavassa imprints a pronounced granulation pattern on the stellar surface. 
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In particular, RSGs give rise to large granules comparable to the 
stellar radius in the H and K bands, and an irregular pattern in 
the optical region. 

This paper is the third in the series aimed at exploring the 
convection in RSGs. The main purpose is to extract photocen- 
tric and photometric predictions that will be used to estimate 
the number of RSGs, detectable by Gaia, for which the paral- 
lax measurement will be affected by the displacements of their 
photometric centroid (hereafter "photocenter"). 

2. RHD simulations of red supergiant stars 

The numerical simulation us ed in th is work has been com- 
puted using CO^BOLD (Frev tag etaD [2QQ2: Freytag & Hofner 
I2OO8I) . The model, deeply analyzed in Paper I, has a mass of 
12 Mo, employs an equidistant numerical mesh with 235^ grid 
points with a resolution of 8.6 Rq (or 0.040 AU), a luminos- 
ity average over spherical shells and over time (i.e., over 5 
years) of L = 93 000 ± 1300 L©, an effective temperature of 
Teff = 3490+13 K, a radius of R = 832 ± 0.7 Rq, and a sur- 
face gravity \ogg = -0.337 ± 0.001. The uncertainties are mea- 
sures of the temporal fluctuations. This is our most successful 
RHD simulation so far because it has ste llar parameters closest 
to Betelgeuse (T ew = 3650 K log ^ = O.O. lLevesque et al.Ll2005l 
or logg = -0.3. lHarper et al.Ll2008l) . We stress that the surface 
gravity of Betelgeuse is poorly known, and this is not without 
consequences for the analysis that will be presented in Sect. 16.11 
(see especially Fig.fTSl). 

For the computation of the intensity maps and spectra based 
on snapshots from the RHD simulations, we used the code 
0PTIM3D (see Paper I) which takes into account the Doppler 
shifts caused by the convective motions. The radiative transfer 
is computed in detail using pre-tabulate d extinction coefficients 
per unit mass generated with MARCS (iGustafsson et al.ll2008b 
as a function of t emperature, density and wavelength for the 
solar composition (lAsplund et al .112006b . The tables include the 
same extensive atomic and molecular data as the MARCS mod- 
els. They were constructed with no micro-turbulence broaden- 
ing and the temperature and density distributions are optimized 
to cover the values encountered in the outer layers of the RHD 
simulations. 



3. Predictions 

In this Section we provide a list of predictions from 3D simu- 
lations that are related to the Gaia astrometric and photometric 
measurements. 

3.1. Photocenter variability 

We computed spectra and intensity maps in the Gaia G band for 
the whole simulation time sequence, namely ~5 years with snap- 
shots ^ 23 days apart. The corresponding spectrum is presented 
in Fig. [Hand the images in Fig.[2l 

Paper II showed that the intensity maps in the optical re- 
gion show high-contrast patterns characterized by dark spots 
and bright areas. The brightest areas exhibit an intensity 50 
times brighter than the dark ones with strong changes over some 
weeks. Paper II reported robust interferometric comparisons of 
hydrodynamical simulations with existing observations in the 
optical and H band regions, arguing for the presence of con- 
vective cells of various sizes on the red supergiant Betelgeuse. 
The Gaia G band images (Fig. ^ are comparable to what has 
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Fig. 1. The transmission curve of the Gaia G band white 
light passband (solid black line), the blue ( dotted blue line) 
and red (dashed red lin e) photometric filters (iJordi et al .1120101: 
iJordi & Carrascol l2007l) together with the synthetic spectrum 
computed from the RHD simulation described in the text. 

been found in Paper II. The resulting surface pattern, though re- 
lated to the underlying granulation pattern, is also connected to 
dynamical eff'ects. In fact, the emerging intensity depends on (i) 
the opacity run through the atmosphere (and in red supergiants, 
TiO molecules produce strong absorption at these wavelengths; 
see spectrum in Fig. [T]) and on (ii) the shocks and waves which 
dominate at optical depths smaller than 1 . 

The surface appearance of RSGs in the Gaia G band aff'ects 
strongly the position of the photocenter and cause temporal fluc- 
tuations. The position of the photocenter is given as the intensity- 
weighted mean of the x-y positions of all emitting points tiling 
the visible stellar surface according to: 

i:f=iZ^.i/o;y)*x(u-) 

Zf.iZf.i/(U>y(U-) 

where / (/, j) is the emerging intensity for the grid point (/, j) 
with coordinates x(i, j), y(i, j) of the simulation, and A/^ = 235 is 
the total number of grid points. Fig. [3] shows that the photocenter 
excursion is large, since it goes from 0.005 to 0.3 AU over 5 
years of simulation (the stellar radius is ^ 4 AU, Fig. [2]). The 
temporal average value of the photocenter displacement is (P) = 
{(Pi -h P^)(i/2)) = 0.132 AU, and or = 0.065 AU. 

At this point, it is important to define the characteristic time 
scale of the convective-related surface structures. RHD simu- 
lations show that RSGs are characterized by two characteristic 
time scales: 

(i) the surface of the RSG is covered by a few large convective 
ceUs with a size of about 1.8-2.3 AU (^60% of the steflar 
radius) that evolve on a time-scale of years (see Fig. |4] and 
Paper I). This is visible in the infrared, and particularly in the 
H band where the H" continuous opacity minimum occurs 
and consequently, the continuum-forming region is more ev- 
ident. 

(ii) In the optical region, as in Fig.O short-lived {a few weeks to 
a few months) small-scale (about 0.2-0.5 AU, ^10% of the 
stellar radius) structures appear. They result from the opac- 
ity run and dynamics at optical depths smaller than 1 (i.e., 
further up in the atmosphere with respect to the continuum- 
forming region). 
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Fig. 2. Maps of the linear intensity (the range is [0 - 230000] erg/s/cm^/A) in the Gaia G band. Each panel corresponds to a different 
snapshot of the model described in the text with a step of about 230 days 5 years covered by the simulation). 
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Fig. 3. Photocenter position, computed from the snapshots of 
Fig. [21 in the Gaia-G band filter. The total simulated time is ^5 
years and the snapshots are 23 days apart. The snapshots are 
connected by the line segments. The dashed lines intersect at the 
position of the geometrical center of the images. Note that the 
photocentric shift stays in the first quadrant for most of the 5 yr 
simulation, and reflects the long lifetime of the large convective 
cell best visible in the infrared H band (Fig. (4]). 




Fig. 4. Map of the linear intensity in the IONIC filter (H 
band as described in Paper I). The range is [0 - 3.1 x 
10^] ergcm"^ s"^ A"^ The snapshot corresponds to the top left 
snapshot at ^ = 21.976 yr in Fig. [2l The large convective cell 
visible in this Figure is swamped in smaller-scale photospheric 
structures in the Gaia G band images. 



Both time scales have an eff'ect on the photocenter excur- the value of dp is mostly fixed by the short time scales corre- 
sion during the 5 years covered by the simulation. On one hand, sponding to the small atmospheric structures. On the other hand. 
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the fact that (Px) and (Py) do not average to zero (according to 
Fig.O the photocenter stays most of the time in the same quad- 
rant, due to the presence of the large convective cell which is 
visible in the H band; see Fig. (Hand Paper I) indicates that the 
5 years period covered by the simulation is not yet long enough 
with respect to the characteristic time scale of the large-scale 
(continuum) cells. 

The top panel of Fig. [5] shows the temporal photocenter dis- 
placement over the ^ 5 years of simulation, which is comparable 
to the total length of the Gaia mission. As seen in the Figure, for 
t < 22 yr, the random displacement is small and increases to a 
maximum value of 0.30 AU at ^ ~ 23 yr. 

In relation with the astrometric implications of this photo- 
centre displacement, which will be discussed in Sect. IH it must 
be stressed that neither (P) nor crp (the latter corresponding to 
the time sampling of the photocentric motion with a rather ar- 
bitrary time interval of 23 days) are the relevant quantities; it is 
instead the standard deviation of P sampled as Gaia will do (both 
timewise and direction wise) which turns out to be relevant. This 
quantity is computed below. 

The bottom panel of Fig. \5\ shows that there is no obvious 
correlation between the photocenter vari ability and the e merging 
intensity integrated in the Gaia G band. iLudwigl (l2006h showed 
analytically that this lack of correlation is to be expected. 

Gaia will scan the sky, observing each object on average 70- 
80 times. The main information that will be used to determine 
the astrometric characteristics of each stars will be the along- 
scan (AL) measurement. This is basically the projection of the 
star position along the scanning direction of the satellite with re- 
spect to a known reference point. By fitting those data through a 
least square minimization, the position, parallax and proper mo- 
tion of the star can be derived. The possibility of extracting these 
parameters is ensured by Gala's complex scanning lavvQ which 
guarantees that every star is observed from many diff'erent scan- 
ning angles. 

In presence of surface brightness asymmetries the photocen- 
ter position will no more coincide with the bary center of the star 
and its position will change as the surface pattern changes with 
time. The result of this phenomenon is that the AL measure- 
ments of Gaia will reflect proper motion, parallactic motion (that 
are modeled to obtain the astrometric parameters of the star) and 
photocentric motion of convective origin. The presence of the 
latter will be regarded as a source of additional noise. 

The impact of those photocenter fluctuations on the astro- 
metric quantities will depend on several parameters, some of 
which are the stellar distance and the time sampling (fixed by 
the scanning law) of the photocentric motion displayed on Fig. [21 
To better a ssess this i mpact, we proceeded as follows. The Gaia 
Simulator (iLuri et alJ 2005) was used to derive scanning angles 
and time sampling for stars regularly spaced (one degree apart) 
along the galactic plane where the supergiants are found. We 
computed the photocenter coordinates at the Gaia transit times, 
by linear interpolation of the photocenter positions of the model 
(as provided by Fig. [3]), after subtracting (Px) (=0.055 AU) and 
(^3;) (=0.037 AU; as we will explain below, a constant photocen- 
tric off'set has no astrometric impact on the parallax). We then 
computed their projection on the AL direction, which we denote 
Pe, being the position angle along the scanning direction on 
the sky. This projection Pg relates to the modulus P of the pho- 
tocenter vector plotted in Fig. [5] through the relation 

P0 = P cos(0-Op) with tan Op = [Py - (Py)) / (Px - (P^)) , (3) 
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Fig. 5. Top panel: modulus of the photocenter displacement as 
a function of time. Bottom panel: photocenter displacement as a 

function of the integrated flux in the Gaia G band, ^25oA^ 
normalized by the temporal average integrated flux. 



and similarly, we define 

Pe = P cos(0 - Op) with tan Op = Py/Px- 



(4) 



^ See http://www.rssd.esa.int/index.php ? project=GAIAj fepage=picture 
_of_the_week&pow= 1 3 



The resulting run of the standard deviation of the photocenter 
displacement with time for two representative stars (one located 
at / = 0° with 59 transits and the other located at / = 241° with 
227 transits) is shown on Fig.[6l which reveals that the time sam- 
pling is, as expected, strongly dependent upon the star position 
on the sky. The transits separated by 2 - 6x10"^ yrs correspond 
to the star being observed in succession by the two fields of view 
(separated by 106.5 degrees on the sky) by the satellite spinning 
at a rate of 6 hours per cycle, whereas the longer intervals are 
fixed by the satellite precession rate. 

Finally, we computed the standard deviation of those pro- 
jections, and obtained ap^ values ranging from 0.06 to 0.10 AU 
(FigO, with {\Pe\} ranging from 1x10""^ to 8x10"^ AU. In the re- 
mainder of this paper, we will adopt ap^ = 0.08 AU. This quan- 
tity, which represents about 2.0% of the stellar radius (^^4 AU; 
Sect. 2), is a measure of the mean photocenter noise induced 
by the convective cells in the model, and it is this value which 
needs to be compared with the Gaia or Hipparcos measurement 
uncertainty to evaluate the impact of granulation noise on the 
astrometric parameters. This will be done in Sections |4] and [51 
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Fig. 6. The along-scan photocenter displacement Pq (in AU) 
against time for two different samplings of the photocenter dis- 
placement of Fig.O corresponding to the Gaia scanning law ap- 
plied to stars located along the galactic plane at longitudes of 
0° and 241°, as indicated on the figures. The top panel of each 
pair provides the distribution of time intervals between succes- 
sive measurements. 



We note that crp^ = 0.08 AU is in fact larger than crp = 
0.065 AU, and this can be understood as follows. First, from 
Eq. (|4]) and basic statistical principles, the following relation may 
be easily demonstrated: 



o-^, = 0.5 ((tI -h {Pfl 



(5) 



under the obvious hypothesis of statistical independence be- 
tween the scanning directions 6 and the photocentric positions 
Pjc, Py. With (P) = 0.132 AU, and crp = 0.065 AU obtained in 
Sect. 13.11 the above relation predicts cr^^ = 0.10 AU, in agree- 
ment with the actual predictions based on Eq. (|4]). If one con- 
siders instead Pg from Eq. © (thus projecting the 're-centered' 
photocentric displacement), there is a small reduction of the 
standard deviation according to 



Crl=CTl.-0.5({P,f^{Py}\ 



(6) 
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yielding crp^ = 0.088 AU, in agreement with the detailed calcu- 
lations shown on Fig. [71 
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Fig. 7. The run of crp^ (expressed in AU; solid blue curve and 
left-hand scale) with galactic longitude for stars located along 
the galactic plane, having a number of transit observations given 
by the red dots (and right-hand scale). 



3.2. Photometric variability 

Another aspect of RSG variability can aff'ect Gaia spectropho- 
tometry. The blue and red photometric bands of Fig. [T] produce 
two spectra of the observed source at low spectral resolution 
(R ^ 50). The photometric system has the advantage of cover- 
ing continuously a wide range of wavelengths providing a mul- 
titude of photometric bands, but it has the great disadvantage 
of being extremely hard to calibrate in flux and wavelength. 
The photometric system of Gaia will be used to characterize 
the star's efl^ecti ve temperature, surface gravity and metallicity 
dih event ^ [2008h . The vigorous convective motions and the re- 
sulting surface asymmetries of RSGs cause strong fluctuations 
in the spectra that will afl'ect Gaia spectrophotometric measure- 
ments (Fig. (8]). In the blue photometric range (top panel), the 
fluctuations go up to 0.28 mag and up to 0.15 mag in the red 
filter (bottom panel) over the 5 years of simulation. These val- 
ues are of the same order as the standard deviation of the visible 
magnitude excursion in the last 70 years for Betelgeuse, 0.28 
mag (according to AAVSCEl)- 

The light curve of the simulation in the (blue - red) Gaia 
color index is displayed in Fig.[9l The temporal average value of 
the color index is (blue -red)= 3.39 + 0.06 at one sigma and there 
are some extreme values at, for example, ^ 22.4 yr, ^ ^ 23 yr, 
and t ^ 25.2 yr. 

Therefore, the unc ertainties on [Fe/H], Teff, and logg given 
by iBailer-JonesI (l2010h for stars with G < 15 should be revised 
upwards for RSGs due to temporal fluctuations from convection. 



3.3. Direct imaging and interferometric observables 

The simulation presented in this work has already been tested 
against the observations at difl'erent wavelengths including the 
optical region (Papers I and II). However, it is now also possi- 
ble to compare the predictions in the Gaia G band to CHARA 
interfer ometric ob servations obtained with the new instrument 
VEGA (Mourard et al. 2009) integrated within the CHARA ar- 
ray at Mount Wilson Observatory. For this purpose, we com- 
puted intensity maps in the blue and red bands of Fig. [T] and 



American Association of Variable Star Observers, www.aavso.org 
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Fig. 8. Spectral fluctuations in the blue and red Gaia photometric 
bands (Fig.© for RSGs: the black curve is the average flux over 
^^5 years covered by the simulation, while the grey shade denotes 
the maximum and minimum fluctuations. The spectra have been 
smoo thed to the Gaia spectral resolution (R ^ 50, iTheveninl 
[20Q8h . 



calculated visibility curves for 36 difl'erent position angles with 
a step of 5° following the method explained in Paper I. FigurefTOl 
shows the intensity maps together with the corresponding visi- 
bility curves. The angular visibility fluctuations are larger in the 
blue band (bottom left panel) because there is a larger contribu- 
tion from molecular opacities (mainly TiO) that shade the con- 
tinuum brightness of the star (top left panel): therefore the sur- 
face brightness contrast is higher. However, in both photomet- 
ric bands the signal in the second lobe, at higher frequencies, is 
-0.2 dex larger than the uniform disk (UD) result, which is mea- 
surable with CHARA. 

The approach we suggest to follow in order to check the re- 
liability of the 3D simulation is the following: to search for an- 
gular visibility variations, as a function of wavelength, observ- 
ing with the same telescope configuration covering high spatial 
frequencies and using the Earth's rotation to study 6-7 diff'er- 
ent position angles in one night. The error bar should be kept 
smaller than the predicted fluctuations: ^^40% in the blue band, 
and ^^20% in the red band at the peak of the second visibility 
lobe. 

To investigate the behavior of the local flux fluctuations, clo- 
sure phases shall bring invaluable information on the asymme- 
try of the source. However, the final consistency check will be 
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Fig. 9. (blue - red) Gaia color index, computed in the blue and 
red photometric filters of Fig.[T] as a function of time. 



an image reconstruction to compare directly the granulation size 
and shape, and the intensity contrast, provided by the planned 
second generation recombiner of the VLTI and CHARA optical 
interferometry arrays. The European Extremely Large Telescope 
(E-ELT, planned to be operating in 2018) with a mirror size 5 
times larger than a single VLT Unit Telescope will be capable of 
near IR observations of surface details on RSGs (Fig.fTTI). 

4. Impact of photocentric noise on astrometric 
measurements 

The basic operating mode of astrometric satellites like Hipparcos 
or Gaia is to scan the sky and to obtain along-scaifl positions rjAL, 
as it was already briefly sketched in Sect. 13. II The core astromet- 
ric data analysis then consists in solving a least- squa re problem 
(for th e sake of simplicity, we neglect the AC term) (iLindegrenI 
l20TQh 



mm 



P,a 



transii 



[7]i - 77(p, a; ti)] 



2\ 



(7) 



for the astrometric parameters p, and the set of satellite attitude 
parameters a, given the transit along-scan positions rji at times 
ti, the model predictions //(p, a; ^/), and the formal error cTj^. on 
the along-scan position /// (including centroiding errors and er- 
rors due to imperfect calibration or imperfectly known satellite 
attitude for instance). If rji is afl'ected by some supplementary 
noise coming from the photocentric motion (which is not going 
to be included in cT;^.), then this photocentric noise of variance 
(Tp^ will degrade the goodness-of-fit in a significant manner, pro- 
vided that CTp^ > (Tn. This statement is easily demonstrated from 
Eq. O, by writing t// = fji -h Pg/., with the first term 77/ represent- 
ing the astrometric motion, and the second term representing the 
along-scan photocentric shift: 



'transit f ~ 



X' 



[fji + Pe,-r]{i^,2i\ti)Y 



(8) 



For the sake of simplicity, we will assume in the following 
that (Tt^ is the same for all measurements. The above equation 
may be further simplified in the case where there is no corre- 
lation between the astrometric and photocentric shifts, so that 



^ Across-scan (AC) measurements will be obtained as well by Gaia 
but will have a lower precision. 




Fig. 10. Top panels: maps of the linear intensity (the range is [0 - 230000] erg/s/cm^/A) computed in the blue and red photometric 
filters of Fig.[TJ Bottom panels: visibility curves from the above maps computed for 36 diff'erent position angles 5° apart (grey lines). 
The black curve is the average visibility while the dashed line is a uniform disk of about the same radius as the simulation snapshot. 
The conversion factor to the more customary unit arcsec"^ on the abscissa axis is arcsec"^ = R^^ • d [pc] • 214.9 (see Paper I). 




proper motiorQ). Although this assumption of absence of corre- 
lation turns out not to be satisfied in real cases (we will return 
to this issue in the discussion of Fig. [2T1) . it nevertheless off'ers 
insights into the situation, and we therefore pursue the analytical 
developments by writing 



X 



^ nrans it trans it 

[^, -77(p,a;f,)f + Pi 



— Xo ~^ Mransit ~ 



(9) 



-4 -2 2 A 

X [AU] 

Fig. 11. Snapshot of 3D simulation (upper left of Fig. O con- 
volved to the PSF of a 42 m telescope (for a star at a distance of 
152.4 pc, see solution #2 in TableO like the European Extremely 
Large Telescope. 



where is the chi-square obtained in the absence of photocen- 

tric motion, and we have assumed o-p^ = ( 1 /transit) ZflT '' ^Ik 
since asymptotically {Pq} = 0. It is important to stress here that 
it is indeed the standard deviation of the photocenter displace- 
ment (sampled the same way as the astrometric data have been) 
- rather than its average value - which matters. In the extreme 
case where there is a constant (non-zero) photocenter shift, there 
will obviously be no impact on the astrometric parameters. 



the cross-product term 2y^rr '(^^ ~ liv^^'^h)) is null. This 
absence of correlation only holds if the photocentric shift oc- 
curs on time scales diff'erent from 1 year (no correlation with the 
parallax), and shorter than a few years (no correlation with the 



^ The large subphotospheric convective cells lead to conspicuous 
spots in the infrared bands, which move on time scales of several years 
(Paper I). However, in the optical bands, these large spots are not so 
clearly visible, since they are swamped in smaller- scale photospheric 
structures. 
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The degradation of the fit due to the presence of the photo- 
centric noise may be quantified through the goodness-of-fit pa- 
rameter F2, defined as 



1/2 



(10) 



where v is the number of degrees of freedom of the;^^ variable. 
The above definition c orresponds to the 'c ube-root transforma- 
tion' of the;^^ variable (IStuart & Ordll 19941) . The transformation 
of (y^, y) to F2 eliminates the inconvenience of having the dis- 
tribution depending on the additional variable y, which is not 
the same for the diff'erent stars. F2 follows a normal distribution 
with zero mean and unit standard deviation. The goodness-of- 
fit F2 thus appears to be an efficient way to detect the presence 
of any photocentric noise. It may be compared to its value F2o 
in the absence of photocentric noise by assuming ;^Q/y = 1 and 
transit ~ 1; then Eq. (O writes 



2 2 



1 -h 



transit ^Pg 



2 A 



(11) 



thus leading to 



F2 = F2o + 



9y 



1/2 



1 + 



1/3 



(12) 



In the case of Gaia, the second term of the above equation 
may be evaluated as a function of crpJcTrj by adopting y = 70, 
as represented on Fig. [121 Since F2 follows a normal distribu- 
tion with zero mean and unit standard deviation, the fit degrada- 
tion will become noticeable if F2 increases by 2 or so, implying 
o'Pglo-j^ ^ 0.6. This translates into a condition on the distance: 



d [kpc] 



CTP, [AU] 
0.6 (T„ [mas] ' 



(13) 



The error on the along-scan position rj should not be confused 
with the end-of-mission error on the parallax (cr^), which ul- 
timately results from the combination of transit transits, with 
transit rangi ng from 59 to 1 20 for Gaia, with an average of 
transit = 78 (lLindegrenll2010l) , and from 10 to 75 for Hipparcos 
(Fig. 3.2.4 of Vol. 1 of the Hipparcos and Tycho Catalogues). 
The number of transits depends (mostly) on the ecliptic latitude. 

For Hipparcos, the individual (Tj^ val ues for each transit may 
be found in the Astrometr ic Data files (Ivan Leeuwen & EvansI 
119981: Ivan Leeuwenll2007ah , and are of the order of 1.7 mas for 
the brightest stars (see Sect.[5]and Fig.O. For Gaia, the quantity 
(Trj may be obtained from the relation 



1/2 



(14) 



where m = 1.2 denotes an overall end-of-mission contingency 
margin, and gpar = 1.91 is a dimensionless geometrical factor 
depending on the scanning law, and accounting for the variation 
of transit across the sky, si nce cr^ is an effective sky-average 
value (see Ide Bruijn32005l). A current estimate of cr^ is 7.8 yuas 
for the brightest stars (lLindegrenl2010l) . yielding cFr^ of the order 
of 30 yuas. To avoid saturation on objects brighter than G - 12.6, 
a special CCD gating strategy will be implemented so that the 
error budget may be assumed to be a constant for G < 12.6 
(Ide Bruiindl2005ULindegrenll2010h . As we show in Sect. [6j only 



I 




Fig. 12. The degradation of the goodness-of-fit F2 in the pres- 
ence of a photocentric motion described in terms of the ratio 
o'Pglo'r^, where cFj^ is the instrumental error. Remember that F2 
follows a normal distribution with zero mean and unit standard 
deviation. 

the bright-star regime matters for our purpose. Inserting these 
values in Eq. ([T3]) , we thus find 



d [kpc] < 0.98 (Tp^ [AU] 
and 

d [kpc] < 55.5 cTp^ [AU] 



for Hipparcos, (15) 



for Gaia. (16) 



Adopting (Tp^ = 0.08 AU for Betelgeuse-like supergiants 
(Sect. l3.ll) yields d < 0.08 kpc for Hipparcos and d < 4.4 kpc for 
Gaia. This limit has to be interpreted as marking the maximum 
distance up to which a photocentric motion with crp^ = 0.08 AU 
will increase the astrometry goodness-of-fit by 2. The validity of 
these conditions will be further evaluated in Sects. [5] and [6] 

In the presence of such photocentric noise, the astrometric 
data reduction process may adopt one of the following three ap- 
proaches: 

(i) Neither the model definition, nor the measurement-error def- 
inition are modified (meaning that the quantities ct/ enter- 
ing Eq. (|9]) are the same as before, and that no attempt is 
made whatsoever to model the granulation). With respect to 
a star with similar properties (same apparent magnitude and 
location on the sky), a star with global- scale convection cells 
will then be recognized by a goodness-of-fit F2 value larger 
than expected depending upon the ratio crpJcTr^ (see Fig.fT^. 
Under those conditions, the resulting /(9rma/ uncertainty on 
the parallax would not be especially large, though; actually, 
it would be exactly identical to the parallax uncertainty in the 
absence of photocentric motion. This is because the formal 
errors on the parameters p (among which the parallax) only 
depend on the measurement errors cFj^ (which were kept the 
same in the presence or absence of photocentric noise), and 
not on the actual measured values (which will be different 
in the two situations). This is demonstrated in Appendix A. 
But of course, the error on the parallax derived in such a way 
is underestimated, as it does not include the extra source of 
noise introduced by the photocentric motion. The next pos- 
sibility alleviates this difficulty. 
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(ii) An estimate of the photocentric noise may be quadratically 
added to cr^ appearing in Eq. ([8]). The error on the par- 
allax will then be correctly estimated (and will be larger 
than the one applying to similar non-convective stars); the 
goodness-of-fit will no longer be unusually large. This is the 
method adopted for the so-called 'stochastic solutions' in the 
Hipparcos reduction, an example of which will be presented 
in Sect. [21 These solutions, called 'DMSA/X', added some 
extra-noise (in the present case: the photocentric noise) on 
the measurements to get an acceptable fit. 

(iii) The model is modified to include the photocentric motion. 
This would be the best solution in principle, as it would al- 
low to alleviate any possible error on the parallax, as they 
may occur with the two solutions above. However, the 3D 
simulations reveal that it is very difficult to model the com- 
plex convective features seen in visible photometric bands 
by a small number of spots with a smooth time behaviour. 
This solution has thus not been attempted. 

The astrometric parameters themselves may change of 
course, for either of the above solution, especially when the pho- 
tocentric motion adds to the parallactic motion a signal having 
a characteristic time scale close to 1 yr. If on the contrary, the 
photocentric motion has a characteristic time scale very differ- 
ent from 1 yr, the photocentric motion averages out, and leaves 
no imprint on the parallax. A similar situation is encountered 
in the presence of an unrecognized orbital motion on top of the 
parallactic motion: only if the unaccounted orbital signal has a 
period close to 1 yr will t he parallax be strongly affected (see 
iPourbaix & JorissenI l200Ql for a discussion of specific cases). 
As explained in Section [TT] RSGs large convective cells evolve 
over time scales of years. In addition, they change slightly their 
position on the stellar surface within the 5 years of simulated 
time (Chiavassa Ph.D. thesi£|), but, unfortunately, it is difficult 
to measure exactly the granule size (Berger et al. 2010, submit- 
ted to A&A) and thus to give a consistent estimation of this dis- 
placement. 



5. A look at Betelgeuse's Hipparcos parallax 

Hipparcos data (lESAl Il997h may hold signatures of global- 
scale granulation in supergiants. The three nearby supergiants 
a Sco (Antares; HIP 80763), a On (Betelgeuse; HIP 27989) and 
a He r (Rasalgethi; HIP 84 345) are ideal targets for this purpose, 
smce lTuthilletal.l(ll997h indeed found surface features on all 
three stars, implying photocentric displacements of the order of 
1 mas (estimated from the product of the fraction of ffux be- 
longing to a bright spot with its radial distance from the geo- 
metric center; see Table [TJ. By chance, observations of the disc 
of Betelgeuse at the time of the Hipparcos mission were done 
by Wilson et al. ( 1992) and Ifuthill et al. (1997) and are shown 
in Fig. [131 iTuthill et al.l reveal that the two bright spots present 
in January 1991 turned into one a year later (January 1992), 
with a much fainter spot appearing at the edge of the extended 
disc. Since the simulation used in this work shows excellent fits 
to the visibility curves, closure phases, and reconstructed im- 
ages based on WHT data in the same filters as those used in 
Fig- El and the fact that RSGs are slow rotators, it is most likely 
that the spots in Fig. [131 are due to convection. Their properties 
have been summarised in Table [H along with the corresponding 
photocentric displacement. The observed photocentric displace- 
ments agree with the model predictions, as can be evaluated from 

^ http://teLarchives-ouvertes.fr/docs/00/29/10/74/PDF/Chiavassa_ 
PhD.pdf 
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Fig. 13. Top panel: interferometric observations of the surface 
of Betelgeuse obtained with the William H erschel Telescope at 
a wavelength of 7 10 nm, in January 1991 (I Wilson et al.lll992h 
and January 1992 (iTuthill et al.lll997h . North is o n top and East 
is to the left. The two images have been taken from lFreytag et al.l 
(2002). The spot properties are summarised in Table [Tl Bottom 
panel: Hipparcos epoch photometry of Betelgeuse with the ver- 
tical lines indicating the epoch of the two interferometric obser- 
vations. 



CTpJmas] = CTpJAU] x ii7[mas] = 0.5 mas, since crp^ = 0.08 AU 
for a supergiant (Sect.O and tu = 6.56+0.83 mas for Betelgeuse 
(see Table O in the remainder of this section, all q uantitie s from 
Hipparcos refers to van Leeuwen's new reduction, [2007bh . 

Considering the fact that the instrumental uncertainty on an 
individual measurement is 1.9 mas for a star like Vega (V = 
0.03) which is as bright as Betelgeuse (V = 0.4 2) (accordin. 
to Veg a's Intermediate Astrometric Data file in Ivan LeeuwS 
'2007a'), Betelgeuse 's convective noise with o-p^ = 0.5 mas 
should be just noticeable on top of the instrumental noise, and 
possibly have some de tectable impact on the astrometric data of 
Betelgeuse. For Vega, Ivan LeeuwenI (l2007al) found a very good 
astrometric solution whose residuals At] have a standard devi- 
ation (Ta// of 1.78 mas, fully consistent with the formal errors 
on T] (top panel of Fig. \T5^. The extreme brightness of Vega 
thus did not prevent from finding a good astrometric solution. 
On the other hand, neither the original Hipparcos processing 
nor van Leeuwen's revised processing could find an acceptable 
fit to Betelgeuse and Antares astrometric data, and a so-called 
'stochastic solution' (DMSA/X) had to be adopted (the kind 
of solution labelled (ii) in the discussion of Sect. |4l), meaning 
that some supplementary noise (called 'cosmic noise') had to be 
added to yield acceptable goodness-of-fit values F2. 

The cosmic noise amounts to 2.4 and 3.6 mas for Betelgeuse 
and Antares, respectively, in van Leeuwen's reprocessing. These 
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Table 1. Properties of the spots observed at the surface of Betelgeuse during the Hipparcos mi ssion. The po s ition offset refers to 
the cen ter of the extended disc, of radius 27 mas (Jan. 1991) and 23 mas (Jan. 1992). Data from I Wilson et al.l (ll992h : lTuthill et all 
([T997h . 
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values correspond to the size of the error bars displayed on 
Figs. \T5\ and [161 Rasalgethi was not flagged as DiVlSA/X, but 
rather as DIMSA/C (indicating the presence of a close compan- 
ion), but its large goodness-of-fit value F2 = 46.63 is indica- 
tive as well of increased noise. Consequently, all three super- 
giants have a parallax standard error larger than expectecQ given 
its Hipparcos magnitude Hp, as revealed by Fig. [14] which dis- 
plays cr^ against the median magnitude for all supergiants (lu- 
minosity classes I and II, of all spectral types) in the Hipparcos 
catalogue. The chromaticity correction has been a serious con- 
cern for the reduct ion of the Hipparcos data of very red stars (see 
iPlatais et al.ll2003l for a detailed discussion of this problem), and 
one may wonder whether the increased noise of the three super- 
giants under consideration could perhaps be related to this effect. 
Since the very red supergiants (with V - / > 2.0) show no appre- 
ciable offset from the rest of the sample in Fig. [14] (at least for 
the brightest supergiants, down to Hp = 8), this possibility may 
be discarded, and the discrepant behavior of Rasalgethi, Antares 
and Betelgeuse in Fig. [14] seems instead related to their large 
apparent brightness, due to their proximity to the Sun. 

Could the poor accuracy of Betelgeuse 's parallax and 
its cosmic noise be related to its s urface features , as al - 
ready sugge s ted in general terms by ^ Barthes & Luri' (1999); 
Gray] (l2000l): iPlatais et al.1 (l2003h : ISvensson & Ludwia (.2005.): 



Bastian & Hefeld (l2005h : lLudwigl([2Q06h : lEriksson & LindegrenI 
(l2QQ7l) . 

The bottom panel of Fig. [15] shows the along- scan residu- 
als Arj for Betelgeuse against time (and Fig. [16] does the same 
for Antares and Rasalgethi), compared with the photocenter dis- 
placements Px and Py determined from the 3D simulation of 
Sect. [HI 

From this comparison, we conclude that the photocentric 
noise, as predicted by the 3D simulations, does account for a 
substantial part of the 'cosmic noise', but not for all of it. A 
possibility to reconcile predictions and observations could come 
from an increase of Betelgeuse's parallax (because the observed 
photocentric motion would then be larger for the crp^ value fixed 
by the models), but this suggestion is not borne out by the re- 
cent attempt t o improve upon Be telgeuse's parallax in the re- 
cent literature (iHarper et al.ll2008l) (solution #2 in Table [2]), by 
combining the Hipparcos astrometric data with VLA positions, 
as this new value is smaller than both the original Hipparcos 
and van Leeuwen's values. The remaining possibility is that the 



^ This larger parallax standard error does not contradict Appendix A 
stating that, in the presence of a photocentric noise, the standard er- 
ror on the parallax should stay the same. This is because this parallax 
standard error is obtained in the framework of a DMSA/X ("stochas- 
tic") solution, where the measurement errors have been artificially in- 
creased by a "cosmic noise" to get an acceptable goodness-of-fit value. 
Hence, the "design matrix" defined in Appendix A, and directly re- 
lated to the variance-covariance matrix of the astrometric parameters, 
has been changed to produce the stochastic solution, thus resulting in 
a larger parallax error. This corresponds to a solution of kind (ii) in 
Sect.a 




10 8 6 4 2 

Median Hipparcos magnilyde 

Fig. 14. Parallax standard errors for supergiants in the Hipparcos 
Catalogue, after van Leeuwen's reprocessing. The dashed line 
represents 1.5 times the standard parallax error for stars with 
precise parallaxes (< 1 0%; solid line) in van Leeuwen's repro- 
cessing (see Fig. 2.19 of I van Leeuwenll2007al) , whereas the solid 
line represents the fiducial relation between the Hipparcos mag- 
nitude and the standard parallax error. Large red circles corre- 
spond to stars redder than V - I = 2. 



3D model discussed in Sect. 13 . 1 l underestimates the photocentric 
motion. In fact. Paper II showed that the RHD simulation fails 
to reproduce the TiO molecular band strengths in the optical re- 
gion (see spectrum in Fig. [TJ. This is due to the fact that the 
RHD simulations are constrained by execution time and there- 
fore use a grey approximation for the radiative transfer. This is 
well justified in the stellar interior, but is a crude approximation 
in the optically thin layers. As a consequence, the thermal gra- 
dient is too shal low and weakens the contrast between strong 
and weak lines (IChiavassa et al.l [20061) . The resulting intensity 
maps look sharper than observations (see Paper II) and thus also 
the photocenter displacement should be aff'ected. As described in 
Paper II, a new generation of non-grey opacities (five wavelength 
bins employed to describe the wavelength dependence of radia- 
tion fields) simulation is under development. This will change 
the mean temperature structure and the temperature fluctuations, 
especially in the outer layers where TiO absorption occurs. 



6. Application to Gaia 

6.1. Number of supergiant stars with detectabie photocentric 
motion 

In this section, we will use Eq. ([T6b to estimate the number of 
supergiants which will have a poor goodness-of-fit as a con- 
sequence of their photocentric motion. This equation requires 
knowledge of crp^, which will be kept as a free parameter in this 
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Table 2. Parallaxes for Betelgeuse computed from various data sets. 
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# Data I^t tu Distance 
[yr] [mas] [pc^ 

1 Hipparcos (FAST+NDAC reduction consortia) 38 2.1 7.63 + 1.64 131.1+^^;^ 

2 van Leeuwen r2007b) 38 2.1 6.56 + 0.83 152.4!^^ ^ 

3 Harper et al. (2008) 51 22.4 5.07 + 1.10 197.2:g-^ 
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Fig. 15. Abscissa residuals A// (in mas on the sky; red dots) al ong 
with the corresponding error bar from Ivan LeeuwenI (l2007ah for 
Betelgeuse and Vega, as a function of time, expressed in years 
from 1991.25. Filled star symbols correspond to the along-scan 
projections Pq of the synthetic photocenter displacements of 
Fig. O The red dashed lines and black solid lines depict the 
±1(T interval around the mean for the Hipparcos data points and 
model predictions, respectively. Note that these displacements 
were computed in the Gaia G filter instead of the Hipparcos Hp 
filter. A test on a given snapshot has shown that the diff'erence 
is negligible: Px = 0.11 AU with the Hp filter, as compared to 
0.13 AU with the G filter. 



section. In Sect. 13.11 crp^ = 0.08 AU was considered as typical 
for Betelgeuse-like supergiants, but Sect. [5] has provided hints 
that 3D models with grey opacities could s omewha t under esti- 
mate thi s quan tity. Moreover, according to iFreytagI (1200 lb and 
iLudwigl (l2006h . crp^ is expected to vary with the star's atmo- 
spheric pressure scale height, which in turn depends upon the 
star's absolute magnitude Mq. To explore the parameter space, 
we thus need to know how crp^ varies with Mq. This is espe- 
cially important since on top of the condition in Eq. ([T6l) relating 
d to o-p^iMc), there is another constraint coming from the re- 
quirement not to saturate the CCD, namely the Gaia magnitude 
G should be fainter than 5.6. All these constraints may be con- 
veniently encapsulated in boundaries in the d - Mq plane, as 
displayed in Fig. [191 

But first, we have to clarify the relation between crp^ and Mq 
which appears to be a critical ingredient in the present discus- 
sion. Unfortunately, 3D hydrodynamical models in the literature 
are scarce. Their main properties are collected in Table O These 
simulations are of two kinds: (i) box-in-a-star modoX^ cover only 
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Fig. 16. Same as Fig.[T5]for Antares and a Her. 



a small section of the surface layers of the deep convection zone, 
and the numerical box includes some fixed number of convective 
cells, large enough to not constrain the cells by the horizontal 
(cyclic) boundaries; (ii) star-in-a-box models, like the one de- 
scribed in this paper (Sect. 2), cover the whole convective en- 
velope of the sta r and have been used to model RSG and AGB 
stars so far (see iFreytag & Hofneill2008l for an AGB model), 
whereas the former simulations cover a large number of stel- 
lar parameters from white dwarf to red giant stars. The transi- 
tion where the box-in-a-star models become inadequate occurs 
around log g ~ 1 , when the influence of sphericity becomes im- 
portant; the star-in-a-box global models are then needed, but 
those are highly computer- time demanding and difficult to run 
so there are only v ery few models available so far. 

ILudwigl (I2006h found that there is a tight correlation between 
the amplitude of the photocentric motion and the size of the 
granular cells. This si ze is related to the pressure scale height 
at optical-depth unity (iFreyt ag 2001.) . The pressure scale height 
is defined as 



mg 



(17) 



where g is the surface gravity, ks is the Boltzmann constant and 
m is the mean molecular mass (m = 1.31 x mn = 1.31 x 1.67 x 
10"24 grams, for temperatures lower than 10 000 K). In the above 
expression, *Hp has the dimension of length. But in the remainder 
of this paper, we adopt instead the simplified definition: 



Teg 



(18) 
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The law relating the standard deviation of the photocenter 
displacement to may be inferred from Fig. [TTl which displays 
the values from Table [S] The transition from the most evolved 
box-in-a-star model (with logH^ ^ 2.57) to our star-in-a-box 
model (log//p ^ 3.85) is still unexplored; consequently, there is 
no guarantee that the trend obtained at log Hp < 2.57 may be ex- 
trapolated to larger Hp values. Different trends are therefore con- 
sidered in Fig. [T71 with a zoom in Fig.[T8l The linear fit of log crp^ 
as a function of log H^ consider s only the box-in-a-star mod- 
els of ISvensson & Ludwigl (l2005l) : the parabolic function is the 
best fit to all the models (including the star-in-a-box supergiant 
model). However, there is strong evidence in the simulations 
that the convective pattern changes strongly from the giant (big 
black circle symbol in Fig. [18]) to the RSG simulations (big black 
squared symbol). The convective related surface structures grow 
enormously in the RSGs and together with the low eff'ective tem- 
perature (i.e., the molecular absorption, strongly related to the 
temperature inhomogenities, is more important) increase the dis- 
placement of the photocenter position (i.e., crp^ is larger). Thus, 
the parabolic fit, which considers all the simulations' configura- 
tions together, is not a completely correct approach because of 
the physical changes reported above. Since the transition region 
between the box-in-a-star (giant stars) and star-in-a-box models 
(RSG stars) is still unexplored, we consider an extreme transition 
by adopting an arbitrary exponential law to relate the last two 
model simulation points (i.e., the transition region between the 
box-in-a-star and star-in-a-box models). Paper I pointed out that 
the reasons for the peculiar convective pattern in RSGs could be: 
(i) in RSGs, most of the downdrafts will not grow fast enough to 
reach any significant depth before they are swept into the exist- 
ing deep and strong downdrafts enhancing the strength of neigh- 
boring downdrafts; (ii) radiative eff'ects and smoothing of small 
fluctuations could matter; (iii) sphericity eff'ects and/or numeri- 
cal resolution (or lack of it). 

To see which among these three possible trends has to be 
preferred, we have made a compilation of photocentric displace- 
ments P from interferometric observations of various super- 
giants available in the literature (see Fig. [18]). Supergiants and 
Miras have been observed several times in the last decade with 
interferometers, often revealing the presence of surface bright- 
ness asymmetries. In several cases (a Ori, a Her, and o Get; see 
Table [4] for the data list and references; more stars will be pre- 
sented in Sacuto et al., in preparation), the observations could be 
represented by parametric models consisting of a uniform disk 
plus one (or more) bright or dark spots. Using the parameters of 
the spots fitting the interferometric data, we computed the posi- 
tions of the photocenter for all observations of a given star and 
from there the standard deviation of these photocentric positions, 
which was then plotted against Hp in Fig. [THJ These observa- 
tional data suggest that the exponential and quadratic fits of the 
simulation data are to be preferred over the linear extrapolation 
of the box-in-a-star values (Fig. [18]). We stress, however, that the 
surface gravity for supergiants like a Her and a Ori are quite 
uncertain (see Table [4]) and also the highly uncertain metallicity 
diff'erences might play a role here. 

The number of stars with photocentric motions detectable 
by Gaia as having bad fits (i.e., large goodness-of-fit F2 values) 
may now be estini ated as follows. The Besan^on Galaxy model 
(iRobin et al.ll2004l) has been used to generate a sample of bright 
giants and supergiants (My < 0) in the region < / < 180 and 
-20 < b < 20 of our Galaxy (where / and b are the galactic 
coordinates). The redd ening has been added se parately using the 
extinction model from iDrimmel et al.l (l2003h . For each one of 
the the 361 069 stars in that sample, we assign the correspond- 




-4-3-2-101234 
log(Hp[s^ K/cm]) 



Fig. 17. Fits to the standard deviation crp^ of the photocentric 
motion predicted from 3D simulations, as listed in Table [3] 
against the pressure scale height Hp. The equation of the (red) 
solid line is logo-p, = -6.110 -h 1.110 log(ifp) with;^^ ^ q.IT: 
th e fit considers only the box -in-a-star models (filled circles) 
of ISvensson & Ludwigl (l2005l) . The (green) dashed line (with 
equation logcr^, = -6.275 -h 1.174 log //p -h 0.039(logifp)2) 
is a fit to all the models of Table [3] (i.e., box-in-a-star mod- 
els and star-in-a-box, the latter being represented by a filled 
square). The (blue) dotted line is an arbitrary exponential law 
that connects the last two points with the following equation 
logcTp^ = -1.09 - 3.434 exp(-0.00149ifp). 
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Fig. 18. Photocenter motions determined from interferometric 
observations for some evolved stars (see Table [4]) overplotted 
on the diff'erent fits of the standard deviation crp^ of the pho- 
tocentric motion as a function of the pressure scale height Hp. 
The large open inverted triangles correspond to the standard de- 
viations of the photocentre deviations for a given observed star. 
Star-in-a-box and box-in-a-star models correspond respectively 
to the large filled square and circle. 



ing expected standard deviation of the photocenter displacement 
(jp^ taken from the exponential or parabolic laws of Fig. [TTl (each 
of these two possibilities being tested separately), with Hp esti- 
mated from Eq. ([T8l)FI 



^ We note in passing that, in the Besan9on sample, there is no star 
matching Betelgeuse parameters if one adopts logg = -0.3 for its sur- 
face gravity, yielding log Hp = 3.85. If on the other hand, one adopts 
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Fig. 19. The location in the (Mg, d) plane for supergiants with a photocenter noise (of standard deviation crp^) significantly altering 
the goodness-of-fit of the astrometric solution (see text), for the case of an exponential link between box-in-a-star and star-in-a-box 
models (see Fig.[T8]and text). 



We then compute the number of stars which fulfill the condi- 
tion expressed by Eq. (O, and having at the same time G > 5.6 
in order not to saturate Gaia CCD detectors. The conversion be- 
tween V and G magnitudes has been done from the color equa- 
tion (adopted from the Gaia Science Performance documen rPh : 

G = y - 0.0107 - 0.0879 (y -/)- 0.1630 (y- 7)^ 

-hO.0086 (y - 7)1 (19) 

With the exponential law, we found 215 supergiants (among 
the 361069 of the full sample, representing half the galactic 
plane) fulfilling these two conditions. They are displayed in 
Fig.[T9]in the d - Mq plane, and are basically confined to a cres- 
cent delineated by the conditions G > 5.6 (corresponding to the 
two lines with an upward concavity, labelled G > 5.6; the two 
lines correspond to two values of the extinction in the G band: 
Ag = and 1) and d < 55.5 crp^iMc) (Eq. ([T6l) , corresponding 
to the green dashed line with a downward concavity). The latter 
line is based on a fiducial relationship between crp^ and Mg, as 

logg = 0.0, we get log Hp = 3.55 and Betelgeuse is then matched 
by stars from the Besan^on sample. This can be seen from the lower 
panel of Fig. [20l since Betelgeuse has Mg = -6.4, when adopting 
^boi - -7.5 from the apparent bolometric flux 1 11.67 x 10~^^ Wcm~ ^ 
dPerrin et al] [2004) and the para llax 6.56 m as jvan Leeuwenll2007ah , 
y - G = 0.98 from y - 7 = 2.32 (lESA|[T997h an d Eq. \l9l BCy =^ -2 05 
from the apparent bolometric flux and V = 0.42 dJohnson et al.lll96^ . 

^ http://www.rssd.esa.int/index.php ?project=GAIA 
&page=Science_Performance 

^ http://www.rssd.esa.int/SYS/docs/lLtransfers/ 
project=PUBDB&id=448635.pdf 



shown on Fig.|20l Some supergiants nevertheless fall outside the 
crescent defined above, simply because of the scatter aff'ecting 
the CTp^ - Mg relationship (Fig. [20|. Obviously, all the super- 
giants of interest are bright in the G band, in the range 5.6 to 
about 8 and will thus be easily identifiable during Gaia data pro- 
cessing. 

With the parabolic law, only one supergiant matches the con- 
ditions: it is the brightest supergiant located in the upper left 
corner of Fig. [20l (green point in the lower panel; note that, in 
Fig. [191 this star is not located below the parabolic threshold line 
as expected, because that line is based on a mean crp - Mq re- 
lation - see Fig. [20]-, and that supergiant happens to have a crp 
value much above average, as seen on Fig. [20|). Thus, Fig. [19] 
suggests that the 'parabolic' link between box-in-a-star and star- 
in-a-box models of Fig. [T7] and [18] is a limiting case: for photo- 
centric motions to be detected by Gaia, the crp^ vs Mq relation 
has to lie above this limiting case (depicted as the green solid 
line in Fig.[T8l). 

In Fig. [19] there is a cluster of stars at Mq = -4.5 (corre- 
sponding to logTeff ~ 3.5 and logg ~ 0.4) which corresponds 
to bright giants or asymptotic giant branch (AGB) stars. They 
are also clearly seen in Fig. [20] as the cluster at crp^ = 0.035 AU 
(with the exponential law) or 0.01 AU (with the parabolic law). 
Since these stars belong to a population diff'erent from super- 
giants (with masses of the order 1 M©), they are not necessar- 
ily confined to the galactic plane as supergiants are. Hence an- 
other sample, now covering a quarter of the sky (0° < / < 180°, 
b > 0°), has been generated from the Besangon model and con- 
tains 702211 giants and bright giants. In this sample, 938 stars 
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Fig. 20. Upper panel: the relation between crp^ and Mg for 
supergiants and bright giants, assuming either an exponential 
(red plusses) or a parabola (blue crosses) to connect the box- 
in-a-star with the star-in-a-box models (Fig. [18]). The equation 
of the dashed green line (through the red plusses) is crp^ = 
-0.29 - 0.10 Mg - 0.0062 , whereas the equation of the 
magenta dotted line (through the blue crosses) is crp^ = 0.083 + 
0.034 Mg + 0.0037 M^. Lower panel: Same as the upper panel, 
for the relation between the pressure scale height Hp and the 
absolute magnitude in the Gaia G band. Only stars with a de- 
tectable photocentric motion (for the exponential fit: red crosses; 
for the parabolic fit: green dot) have been plotted. 



satisfy the condition of detection of the photocentric motion with 
the exponential law, and none with the parabolic law. The rela- 
tion Mg-ctpq thus appears as an essential ingredient, but unfortu- 
nately quite uncertain still, especially for those among the bright 
giants which are pulsating as long-period variables. The pulsa- 
tion makes the niodelling especially difficult (se e for instance 
iFrevtag & Hofned 120081: IChiavassa et"aD l2010bl for an appli- 
cation of 3D AGB models to the star VX Sgr). Nevertheless, 
numerous observations have revealed their surface brightness 
asymmetries (e.g.. iRagland et al.ii2006i. and references therein). 



6.2. Impact on the parallaxes 

To evaluate the impact of the photocentric shift on the parallax, 
we proceeded as follows. The sampling times, scanning angles. 



along- scan measurements and their err ors were obtaine d from 
the Gaia Object Generator v7.0 rGOGR: llsasi et al.l2010h for the 
supergiant stars from the sample generated using the Besan9on 
model described in the previous section. A photocentric motion 
deduced from the photocentre position computed from the snap- 
shots of the red supergiant model (see Fig. [3]) was added on the 
along scan measurements (the photocentric shift was converted 
from linear to angular shifts, according to the known stellar dis- 
tance). The red supergiant model gives a single photocentre po- 
sition sequence. Yet the sequence for every star should be diff'er- 
ent. Therefore the sequence was rotated for every star by a ran- 
dom angle before being added to the along scan measurements. 
The astrometric parameters were then retrieved by solving the 
least-squares equation (Eq. d?])), separately with and without sur- 
face brightness asymmetries. The resulting parallaxes are com- 
4 pared in Fig. [21] 

Fig.[23]presents the histogram of the quantity {tu-tu^^o^IcTjj 
for three diff'erent ranges of distances. It is clearly seen that the 
distribution, quite peaked at zero for distant stars, becomes wider 
for nearer stars, meaning that the ratio of the error on the parallax 
to its formal error increases with decreasing distance. Similarly, 
the fits of the astrometric data are worse for stars closer by, 
and this eff'ect is clearly seen on Fig. [24] displaying the rela- 
tion between the goodness-of-fit parameter F2 and the distance. 
The run of F2 with distance is consistent with that predicted by 
Eq. O, for y = 70, ap^ = 0.1 AU and CTr^ = 0.03 mas. 

Coming back to Fig. [22] it is remarkable that the relative er- 
ror on the parallax, namely (ru - ^spot)/^ is almost independent 
of the distance and amounts to a few percents. This is in fact 
easy to understand, if one assumes that the diff'erence tu - Tu^^oi 
must somehow be proportional to the amplitude of the excur- 
sion of the photocenter on the sky, which must in turn be re- 
^lated to 6, the angular radius of the star on the sky; therefore, 
{tu - ^spot)/^ = c^OItu = aR, where a is the proportionality 
constant and R is the linear radius of the star (expressed in AU). 
Thus we conclude that the relative error on the parallax is inde- 
pendent of the distance, and is simply related to the excursion of 
the photocenter expressed in AU. 

These simulations for a sample of Betelgeuse-like super- 
giants thus allow us to confirm the results obtained in Sect. 16.11 
(and Fig. [19]), in particular the fact that the impact on the 
goodness-of-fit remains noticeable up to about 5 or 6 kpc 
(Fig. [24]). 



7. Conclusion 

We have provided astrometric and photometric predictions from 
3D simulations of RSGs to evaluate the impact of the surface 
brightness variations on the astrometric parameters of these stars 
to be derived by Gaia. 

We found that the global- scale convective pattern of RSGs 
cause strong variability in the position of the photocenter, P. 
From a 3D simulation of a Betelgeuse-like supergiant, (P) = 
0.132 + 0.065 AU (i.e., more than 3% of the stellar radius) show- 
ing excursions from 0.005 to 0.3 AU over the 5 years of sim- 
ulation. In addition, the spectra show large fluctuations in the 
red and blue Gaia bands of up to 0.28 mag in the blue and 0.15 
mag in the red. The Gaia color index (blue - red) also fluctu- 
ates strongly with respect to time. Therefore, the uncertainties 
on [Fe/H], Tq^ and logg should be revised upwards for RSGs 
due to their convective motions. We have furthermore provided 
predictions for interferometric observables in the Gaia filters that 
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Table 3. Photocenter motion from 3D simulations. is given with two different dimensions. 



Model 


Configuration 


log//p 


logH 


log(crp^) 


^efF 


logg 


R 


Mboi 






(s^K/cm) 


(10^ cm) 


(AU) 


(K) 




(Rg) 




White dwarfs 


box-in- a- star 


-3.92 


-3.15 


-10.28 


12000 


8.00 


1.28 10-' 


11.03 


Sun^ 


box-in- a- star 


-0.68 


0.12 


-7.12 


5780 


4.44 


1 


4.74 


Procyon A"^ 


box-in- a- star 


-0.19 


0.61 


-6.37 


6540 


4.00 


2.10 


2.59 


^ Hydrae^ 


box-in- a- star 


0.75 


1.55 


-5.45 


4880 


2.94 


10.55 


0.36 


Cepheid^ 


box-in- a- star 


1.66 


2.46 


-4.25 


4560 


2.00 


30.17 


-1.63 


Red giant"^ 


box-in- a- star 


2.57 


3.36 


-3.02 


3680 


1.00 


95.25 


-3.19 


RSG^ 


star-in- a-box 


3.88 


4.68 


-1.10 


3490 


-0.34 


832 


-7.66 



" iSvensson & Ludwid (120051) 



^ This work, Sect.O 



Table 4. References used to compute the photocentric shifts from interferometric data. Stellar parame ters are from lLevesque et aP 
(l2QQ5h a ndlHarper etal l (l2008h for a Ori and a Her (set 2). lEl Eidl(ll994l) for a Her (set 1) and lFernid (Il995h for o Get. Parallaxes 
are from lvan LeeuwenI (l20Q7bl) . Only observations in the optical range have been kept. 



Name log g 



(K) 



log//p 
log (s^ K/cm) 



(mas) 



(mas) (AU) 



A 

(nm) 



Date References 



I Harper et al.l (I2008|) f s eFT) 
LevesQue etal. ( 20051) (s et 2) 
Busche r et al l I199qI), IWilson 
Tuthill et all 119971). lYoung et al 
(120071) . iHaubois et all (1200^ 



a Ori 



-0.3 
0.0 



3650 
3650 



3.86 
3.56 



a Her 



0.76 
0.0 



3400 
3450 



2.77 
3.53 



o Get -0.6 2900 



6.56 
6.56 



9.07 
9.07 



4.06 10.91 



1.216 0.185 700 02/1989 



1.637 


0.249 


710 


01/1991 


0.369 


0.056 


700 


01/1992 


0.694 


0.106 


700 


01/1993 


0.550 


0.084 


700 


09/1993 


0.427 


0.065 


700 


12/1993 


0.144 


0.022 


700 


11/1994 


0.395 


0.060 


700 


12/1994 


0.302 


0.046 


700 


12/1994 


0.142 


0.021 


700 


01/1995 


0.025 


0.004 


700 


01/1995 


0.009 


0.001 


700 


11/1997 


(P) 


0.075 








0.011 






o-p 


0.075 






0.340 


0.037 


710 


07/1992 


0.765 


0.084 


710 


06/1993 


(P) 


0.060 






iP') 


0.004 






o-p 


0.033 






1.202 


0.110 


710 


07/1992 


0.850 


0.078 


700 


01/1993 


0.990 


0.091 


710 


09/1993 


1.950 


0.179 


710 


12/1993 


(P) 


0.114 






(P') 


0.015 






o-p 


0.045 







^TZI ([19921. [T997I) 
' (I200d) . iTatebe et all 



El Eid (1994) (set 1) 
Levesque et al. (20 05) (set 2) 
lTuthilletal.l ( [T^ 



iTuthill et all ( [199^ 



can be tested against observations with interferometers such as 
VEGA at GHARA. 

Then we studied the impact of the photocentric noise on the 
astrometric parameters. For this purpose, we considered the stan- 
dard deviation of the photocenter displacement predicted by the 
RHD simulation, sampled as Gaia will do (both timewise and 



direction wise). We called this quantity crp^, where 6 is the po- 
sition angle of the scanning direction on the sky, and we found 
CTp^ = 0.08 AU for Betelgeuse-like supergiants. This photocen- 
tric noise can be combined with cTr^ = 30 yuas (the error on the 
along-scan position rj) for Gaia to determine the maximum dis- 
tance (d < 4.4 kpc) up to which a photocentric motion with 
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Fig. 21. Comparison of parallaxes for supergiant stars with and 
without surface brightness asymmetries (spots), normalised to Fig. 23. Histograms of the relative error on the parallax of su- 
o-^. The stars falling on the horizontal line with ordinate are pergiant stars, for different ranges of distances, 
very reddened stars, are consequently quite faint, and therefore 
have large errors on their astrometric measurements and thus on 
their parallax. 



B 
B 




3 4 5 

distance [kpc] 



Fig. 22. 

how the 
tance. 



Same as Fig. [2T1 but normalized by the parallax. Note 
relative parallax error is almost independent of the dis- 
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Fig. 24. Same as Fig. [21] for the goodness-of-fit F2. The solid 
line corresponds to the prediction from Eq. ([T2l) , with y = 70, 
cTp^ = 0.1 AU and ctj^ = 0.03 mas. Part of the large scatter at 
intermediate distances (2 to 4 kpc) is due to differ ent CCD gating 
sequence (see Ide Bruiindl2005l: lLindegrenll2010l) . 



CTpg = 0.08 AU will generate an astrometric noise of the order 
of the astrometric error on one measurement (more precisely 0.6 
times that error, yielding an increase of the F2 goodness-of-fit 
parameter by 2 units). The value crp^ = 0.08 AU could even 
be somewhat underestimated, as we guessed from the compar- 
ison of the along-scan Hipparcos residuals for Betelgeuse with 
the RHD predictions. We concluded that the predicted photocen- 
tric noise does account for a substantial part of the Hipparcos 
'cosmic noise' for Betelgeuse and Antares, but not for all of it. 
This may be due to the fact that the temperature stratification in 
the RHD models is not completely correct due to the grey ap- 
proximation used for the radiative transfer. The resulting inten- 
sity maps have higher contrast than the observations, as shown 
in Paper II, and the photocenter position can thus also be af- 
fected. New simulations with wavelength resolution (i.e., non- 
grey opacities) are in progress and they will be tested against 
these observations. 



We estimated how many RSGs might have have an abnor- 
mally large goodness-of-fit parameter F2. We found that the 
photocentric noise should be detected by Gaia for a number 
of bright giants and supergiants varying between 2 and about 
4190 (215 supergiants in each half of the celestial sphere and 
940 bright giants in each quarter of the sphere; see Sect. 16.11) , 
depending upon the run of crp^ with the atmospheric pressure 
scale height Hp, and to a lesser extent, depending on galactic 
extinction. The theoretical predictions of 3D simulations pre- 
sented in this work will be tested against the multi-epoch in- 
terferometric observations of a sample of giants and supergiants 
(Sacuto et al. in preparation), with the hope to better constrain 
this (Tp^ - Hp relation. In a forthcoming paper (Pasquato et al., 
in preparation), we will evaluate how the Gaia reduction pipeline 
behaves when facing the bright-giants and supergiants granula- 
tion. More specifically, we will show that the distance to the star 
is the main driver fixing which one among all the possible solu- 
tion types (single-star, acceleration, orbital, stochastic) is actu- 
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ally delivered by the pipeline (the acceleration and orbital solu- 
tions being obviously spurious). 

Finally, a very important conclusion is that the parallax for 
Betelgeuse-like supergiants may be affected by an error of a few 
percents. For the closest supergiants {d < 2.5 kpc), this error 
may be up to 15 times the formal error cr^ (see Fig. [21]) resulting 
from the measurement errors and estimated from the covariance 
matrix. In a forthcoming paper (Pasquato et al., in preparation), 
we will moreover show that this error is sensitive to the time 
scale of the photocentric motion (which is in turn fixed by the 
granulation and the stellar rotation). 

There is little hope to be able to correct the Gaia parallaxes 
of RSGs from this parallax error, without knowing the run of 
the photocentric shift for each considered star. Nevertheless, it 
might be of interest to monitor the photocentric deviations for 
a few well selected RSGs during the Gaia mission. Ideally, this 
would require imaging the stellar surface, although monitoring 
of the phase closure on three diff'erent base lines may already 
provide valuable information on the size of the inhomogeneities 
present on the stellar surface (see Sacuto et al., in preparation). 
The best suited targets for that purpose would be supergiants 
with G magnitudes just above the Gaia saturation limit of 5.6, 
where the astrometric impact is going to be maximum, and at 
the same time, still within reach of the interferometers. The cor- 
responding diameter will be on the order of 4 mas (derived from 
the radius 8307?© for a Betelgeuse-like supergiant seen at a dis- 
tance of 2 kpc if G = 5.6, Ag = 1, and Mq = -6.6). A search 
for G, K or M supergiants (of luminosity classes I, la, lab or lb) 
with 5.6 < y < 8 in the SIMBAD database yielded only three 
stars (XX Per, HD 17306 and WY Gem) matching these criteria, 
the latter being a spectroscopic binary which will disturb the ra- 
dius measurement and is thus unsuited for this purpose. It may 
therefore be necessary to select such targets from the Gaia data 
themselves, after the first year of the mission. 
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Appendix A: Formal errors on the parameters of a 
least-squares minimisation 

We provide here a shor t demonstration o f a well-known sta- 
tistical result (see e.g., iPress et al]|l992h . which may appear 
counter-intuitive in the present context, namely the fact that the 
presence of an extra- source of unmodelled noise will not change 
the formal errors on the parameters derived from a least- squares 
minimisation. 

Consider the case where the data points (xi,yi) (i = 1, ...N) 
must be fitted by a general linear model 

y(x) = ^^^,akXk(x) (A.l) 



18 



A. Chiavassa et aL: Spectro-photocentric and photometric variability of red supergiant stars 



where Xk(x) (k = I, ...M) are M arbitrary (but known) functions 
of X, which may be wildly non linear. The merit function is de- 
fined as 



X 



E 

/=1 



(A.2) 



where cr/ is the measurement error on yi presumed to be known. 
To simplify the notation, we define the design matrix A (of size 
X M) by 



Xj(xi) 



(A.3) 



the vector b of (normalized) measured values, of length A^: 



yi_ 



(A.4) 



and finally the vector a of length M whose components are 
the parameters ak (k = I, ...M) to be fitted. The least-squares 
problem may thus be rephrased as 



find a that minimizes = |A • a - bp, 
whose solution may be written 

(a'^ • a) • a = A^ • b, (A.5) 

with C = (A^ • A^ being the variance-covariance matrix de- 
scribing the uncertainties of the estimated parameters a. The 
crucial point to note here is the fact that matrix C involves the 
measurement uncertainties ct/ but not the measurements yi them- 
selves. Therefore, changing yi, in the presence of an unmod- 
elled process (like photocentric motion) without changing the 
measurement uncertainties ct/, will not change the formal errors 
on the resulting parameters a. But of course, along with the 
goodness-of-fit parameter F2 (see Eq. (fTOl) ) will be larger in the 
presence of a photocentric noise, as the scatter around the best 
astrometric solution will be larger than expected solely from the 
measurement errors. Therefore, it is F2 and its associated 
but not the formal parallax error, which bear the signature of the 
presence of photometric noise. 



In fact, this statement only holds in the case where the errors ct/ 
are normally distributed, which is supposed to be the case for the spe- 
cific problem under consideration (namely, the Gaia along- scan mea- 
surement errors). 



